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Abstract 

Induced by three gluons symmetry, s, t, u symmetric expressions are widely involved in collider 
physics, especially in heavy quarkonium physics. In this work we study general form of s, t, u 
symmetric polynomials, and find that they can be expressed as polynomials where the symmetry 
is manifest. The general form is then used to simplify expressions which asymptotically reduces 
the length of original expression to one-sixth. Based on the general form, we reproduce the exact 
differential cross section of J ftp hadron production at leading order in v 2 up to four unknown 
constant numbers by simple analysis. Furthermore, we prove that differential cross section at 
higher order in v 2 is proportional to that at leading order. This proof explains the proportion 
relation at next-to-leading order in v 2 found in previous work and generalizes it to all order. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) is currently believed to be the fundamental theory of 
the strong interaction. Thanks to asymptotic freedom and factorization |l|-|6| properties, ap- 
plication of QCD to a physical process can be factorized into convolution of non-perturbative, 
but universal, long distance matrix elements with perturbative calculable infrared-safe 
short distance coefficients. Short distance coefficient at each order in perturbative expansion 
can be typically understood as the differential cross section of two partons scattering to pro- 
duce n partons. This 2 — > n parton level process is conveniently obtained from the process 
— > 2 + n through crossing. Although partons can be either gluon or (anti-)quark, in many 
cases these 2 + n partons are mainly gluons which are the gauge bosons of QCD. Therefore, 
the identity property between gluons results in a large symmetry for these processes. 

Considering the specific type of processes — > g{ki) + g{k,2) + g(k 3 ) +H(P), where k { are 
momentum of each gluons, P is the total momentum of a "cluster" H which includes one or 
more partons. The "cluster" here means differential cross section of the process should be 
sensitive only to the total momentum P of the cluster, but insensitive to the detail within 
the cluster. Momentum conservation gives 

k% + k£ + k£ + P' i = 0, (1) 

with kf = and P 2 = M 2 , where M is the invariant mass of cluster H. The associated 
Mandelstam variables are given by 

s = (h + k 2 ) 2 = (P + k 3 ) 2 , (2a) 
t = (k 2 + k 3 ) 2 = (P + k 1 ) 2 , (2b) 
u = (k 3 + h) 2 = (P + k 2 ) 2 , (2c) 

with s + t + u = M 2 . Symmetry between the three gluons implies that the following function 
is symmetric under the exchange of s, t and u: 

F(M 2 ,s,t,u)=^\A\ 2 , (3) 

where A is Feynman amplitude of the process and ^2 means summation/average over spin 
and color of these three gluons. Before doing any integration, F(M 2 ,s,t,u) is a fraction 
polynomial, and both its numerator and denominator are symmetric polynomials. As a 
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result, to study the restriction introduced by gluons symmetry, it is equivalent to study the 
property of symmetric polynomial. 

The above type of processes are widely involved in collider physics. The simplest example 
is studying two jets production in hadron colliders and calculation of four gluons scattering 
is need, where H is also a gluon. In heavy quarkonium physics, there are a lot of processes 
belong to this type, including heavy quarkonium decays to light hadrons [PHT3]] and heavy 
quarkonium production in hadron colliders jL4j-[l8[], where if is a heavy quark anti-quark 
pair with a very small relative momentum. Hopefully, study the property of symmetric 
polynomial will introduce rigorous restriction for these processes. 

Among others, there is a very interesting finding in heavy quarkonium physics recently 
that, at leading order in a s and in the large transverse momentum p? limit, relativistic 



correction term for J/ip hadron production is proportional to the leading term [17]. The 
proportion behavior is nontrivial because there are more than one parameters even in the 
large transverse momentum limit. It is likely to have a symmetry to protect this behavior, 
which is one of the motivations to study the symmetry property induced by three identical 
gluons. 

The rest of the paper is organized as follows. We study the general form of a s, t, u 
symmetric polynomial in Sec. [TJ. In Sec. 11 A we devote to massless case and find the poly- 
nomial can be expressed in a form where symmetry is manifest. By explicitly constructing, 
we generalize the massless result to include also massive particles in Sec. 11 B . The general 
form is then used to simplify expression in Sec. Asymptotically, this simplification can 
reduce the length of a expression to one-sixth. In Sec. [IV], we use the general form to re- 
produce some known results and explain the unexpected proportion behavior of relativistic 
correction of J/ip hadron production |l7[j- The proportion behavior is also generalized to 
all order in v 2 . Finally, we summarize the results in this work and give a outlook for future 
works in Sec. [V[ 



II. GENERAL FORM OF s, t, u SYMMETRIC POLYNOMIAL 

In this section, we study the general form of a polynomial F n (M 2 , s,t,u) which is sym- 
metric under the exchange of s, t and u. Because all terms in F n (M 2 ,s,t,u) should have 
the same mass dimensions, a explicit subscript n is attached to denote the mass dimensions 
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2n. Define three symmetric combinations 

Si :—s + t + u, (4a) 
S 2 ;=±(s 2 + t 2 + u 2 ), (4b) 
S 3 :=^(s 3 + t 3 + u 3 ). (4c) 

We will find that F n (M 2 , s, t, u) can be expressed as a polynomial F n (Si, S 2 , S3). This result 
is not hard to be understood because the relation s + t + u = M 2 has already suggested 
that only three variables are independent in F n . What is not clear is that whether F n is a 
polynomial. In the following, we first study the case where M — 0. Then the extension to 
massive case is straightforward. Main results in section are Eqs. (fTS]) and (|T§D. 

A. Massless case: s + t + u = 

Define 

f n (s,t,u) := F n (0,s,t,u) , (5) 
which is a symmetric homogeneous polynomial. The general expression of f n (s,t,u) is 

n+l 
2 

f n ( S , t, U ) = J2 X i [^(^ + Ul ) + tn ~\ U% + S ') + un '\ si + t 1 )} + Stu fn-3(s, t, U) , (6) 

8=0 

where Xi are independent of s, t and u, f n s(s,t,u) is a symmetric polynomial with power 
n — 3. Obviously, f n _ 3 (s,t,u) = if n < 3. Note that in this case Si = s + t + u = 
and S3 = |(s 3 + 1 3 + m 3 ) = stu, therefore, we need only prove the following function is a 
polynomial in S2 and S3: 

n+l 
2 

/ n (s, U) = J2 X i [^"'^ + U ') + + ^ + ^"(^ + **)] • ( ? ) 

We prove it recursively. 

• n = or 1 is trivial. 

• n = 2: 

/ 2 (s, t, u) = 2x (s 2 + t 2 + u 2 ) + 2a;! (st + tu + us) = 2(2x - Xt)S 2 . (8) 
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n > 3 and n G 2iV + 1: 

In this case, we will show that f n (s,t,u) oc s, and then using symmetry we can get 
f n (s,t,u) oc stu, that is f n (s,t,u) = S 3 f n ^ 3 (s,t,u). Setting s = 0, one has t = —u, 
thus 



n+l 
2 



n + l 
2 



/ n (o, = [r-v + = + (-1)'] u n = . (9) 

i=0 i=0 

Recall the general expression of f n (s,t,u) in Eq. (0), the above result means 
f n (s,t,u) oc s. 

n > 3 and n G 2N, but n 6N: 

In this case, we will show that f n (s,t,u) oc S2. Solutions of 



is 



51 = 0, 

s 2 = o, 

s = e 3 u 



t = u . 



(10) 



(11) 



Thus, equivalently, we need to show that f n (s,t,u) vanishes for solutions in Eq. ( |TT 

fnie^u, e Ti %u, u) 



--11" > x, \e^" 3 ( e +t 3 

j=0 



[ 

3=0 



■ 2(n-j)-K r _y2j 



l) + e^^3^(e ±l 3 +l) + (e 



3^' 3 _|_ gT' 3 



;2jTV 



e 3 + e 3 _)_ e il 3 + c.c. 



i=o 

2 



(12) 



3=0 



2(n-2j)n 2(n-j> , 2j7T 
cos h cos h cos — — 



Because cos 2 -f = 1 if j G 3N and cos 2f = -§ if j g 3iV, cos 2(n ~ 3 2j> + cos + 
cos = for all j, that is f n (e ±l ^~u, e Tl ^~u, u) = 0. 

n > 3 and n G 6iV: 
Defind 



Gi(s, t, u) := s n -\f + ij) + t n " V + s*) + u n_ V + f ) . 



(13) 



It is easy to find that 



Gj(0,-u,u) = 2(-l)% n , 



G 3 (e 



3 u, e + 3 u, u) 



2jTT 



(14) 



Qe l 3 u 



Thus 



e s G,, (e 3 u, e + s u,u) — e 



. 2 jo 7T 



32 



e 3 e + 3 u, u) 



(-1)^(0, -u,u) - (-1)^(0, -«,«) = 0, 
-G n (e 

which means 

(-l) jl G n (s,t,u) - (-l)^G h (s,t,u) oc S 3 
pr l ^3~Gj 1 (s, t, u) — e~ l ^i~Gj 2 (s, t, u) oc S 2 ■ 



(15) 



(16) 



For any ji, there exists a j 2 which guarantees the coefficient matrix of Eq. ([Tj]) to be 
non-zero. Therefore, solution of Eq. (|16|) gives 



Gj(s, t, u) = S 2 Aj(s, t, u) + S 3 Bj(s, t, u) 



(17) 



where Aj(s,t,u) and Bj(s,t,u) are symmetric polynomial. Specifically, taking advan- 
tage of results in previous cases, we find Aj(s,t,u) oc S 2 and Bj(s,t,u) oc S3. Using 
this argument recursively, one gets that Gj(s, t, u) is a polynomial in S 2 and Sf . As a 
result, f n (s,t,u) is a polynomial in Sf and Sf. 

Combine all possible cases above, we indeed proved that f n (s,t,u) is a polynomial in S 2 
and S3 for any n. More precisely, we find 

X f (S 2 3 ,S 3 2 ),nG6iV, 
SlSsXn^r (Sf, Sl),ne6N + l, 

6 

5 2 I»-2 (S|, S|), n G 6iV + 2 , 

~ 3 (18) 

S 3 A^3 (Sf , S|), n G 6iV + 3 , 

6 

S 2 Xn^4 (Sf , S 3 2 ), n G 6iV + 4 , 

6 

S 2 S 3 Xu^ (S'l Si), n G 6N + 5 , 

6 

where Xi(x,y) is an arbitrary homogeneous polynomial in x and y with power i. Specially, 
Xi(x,y) = if i < 0. 



f n (s,t,u) = f n (S 2 ,S 3 ) = < 



B. Massive case: s + t + u = M 2 

When Si = M 2 ^ 0, we will show that F n (M 2 , s, t, u) can be expressed as: 

n 

F n (M 2 , s, t, u) = F n (S u S 2 , S 3 ) = ^2 S\f n ^{S 2 , S 3 ) , (19) 

i=0 

where general form of f n ~i(S 2 , S 3 ) is explicit in Eq. (|T^). We prove Eq. (|19j) by construction. 

First, i — term in Eq. ( |T9"D can be obtained by setting Si = M 2 = in F n (M 2 , s, t, u) 
using the method discussed in Sec. [11 A| . Then we find F n (M 2 , s, t, u) — f n {S2, S3) is zero if 
one sets Si = 0, that is F n {M 2 , s,t,u) — fni^jSs) oc Si. Considering that F n [M 2 , s,t,u), 
fn(S2, S3) and 5*i are symmetric polynomials, 

Fn _ l( M>, ., t, u) = ^-^-MSi.S,) 

Si 

is also a symmetric polynomial. Applying this method repeatedly, result of Eq. ( |19"D can be 
achieved. 

III. SIMPLIFY EXPRESSIONS 

General forms in the last section can be used to simplify expressions that are symmetric 
polynomials in s, t and u by expressing them in terms of manifest symmetric form. In large 
n limit, one finds from Eq. (|l^) that the symmetric form is a polynomial with power ^, 
therefore, asymptotically this method can reduce the length of original expression to one- 
sixth. Also, the method can be easily realized in terms computer program. Here, we give 
two examples. 

The first example is an ideal expression 

42 

/ 42 : = s\t A2 - 1 + u A2 - 1 ) + tV 2 "' + s A2 ~ l ) + u^s 42 ^ + t A2 - 1 ) , (21) 
i=i 
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with s + t + u = 0. Its explicit expression in terms of t and u is 

/ 42 = 2(t 2 + + w 2 ) 3 [t 36 + 18t 35 w + 342t 34 w 2 + 4029t 33 w 3 + 34542t 32 w 4 + 229536£ 



3 V 



+ 1229611t 30 M 6 + 5455041t 2 V + 20431896t 2 V + 65541187t 2 V + 182032158t 2 V° 
+ 441440337t 2 V 1 + 940899497t 2 V 2 + 1771715799t 2 V 3 + 2959077438t 2 V 4 
+ 4396930001t 21 w 15 + 5825638020t 2 V 6 + 6892901679t 1 V 7 + 7289748245t lf V 8 
+ 6892901679^ V 9 + 5825638020t 1 V° + 4396930001t 1 V 1 + 2959077438t 14 w 22 
+ 1771715799t 13 u 23 + 940899497t 1 V 4 + 441440337t n w 25 + 182032158t 1 V 6 
+ 65541187tV 7 + 20431896tV 8 + 5455041tV 9 + 1229611tV° + 229536t 5 M 31 
+ 34542tV 2 + 4029tV 3 + 342tV 4 + 18tu 35 + u m ] . 

(22) 

As is symmetric under exchanging of s, t and u, we can express it in terms of S2 and S3. 
Specifically, we find result for terms within brackets of Eq. (^) is 

= + mS^S 2 + 3060S 2 12 S 4 + 12376S 2 9 S 3 6 + 12870S ,6 S 3 8 + 3003S 3 S 3 10 + 91S 12 , 

(23) 

which is much simpler than the original expression. Result in Eq. ( PB"D has two meanings. 
Firstly, it tests the proof in Sec. |11 A| to be true by explicit calculation. Secondly, it shows 



the asymptotic behavior that our method can reduce the length of original expression to 
one-sixth. 

The second example is a massive one. We take the equation (A5d) of Ref. |15 
\A(gg -> ccfp' 8 ^)! 2 = ^p-^— L_£l____ j 

+ s 2 /(s 2 - /) 4 + M 2 sz 2 (s 2 - /) 2 (3s 2 - 2/)(2s 4 - 6s 2 / + 3/) 
+ M 4 [9s 12 - 84s 10 / + 265s 8 / - 382s 6 / + 276s 4 / - 88s V° + 9/ 2 ] 

- M 6 s[54s 10 - 357s 8 / + 844s 6 / - 898s 4 / + 439s 2 / - 8l/°] 

+ M 8 [153s 10 - 798s 8 / + 1415s 6 / - 1041s 4 / + 301s 2 / - 18/°] ^ 

- M 10 s [270s 8 - 1089s 6 / + 1365s 4 / - 616s 2 / + 87/] 
+ M 12 [324s 8 - 951s 6 / + 769s 4 / - 189s 2 / + 9/] 

- 9M 14 s(6s 2 - /)(5s 4 - 9s 2 / + 3/)3M 16 s 2 (51s 4 - 59s 2 / + 12/) 
-27M 18 s 3 (2s 2 -/) + 9M 20 s 4 }, 



which is the squared amplitude of 3 Pg channel for J / ip hadron production via gluon-gluon 
fusion. Even a new variable z := ytu is introduced, the expression is still very long. Using 
our method, it can be expressed in a much compact form. Terms within braces of Eq. ( |24| ) 
can be expressed as 

{•••} = 5 2 4 5 2 + SxSlSz (6S| - Si) + Si (95 2 6 + 25 3 5 2 + 55 3 4 ) + 25 3 5 2 5 3 (65 3 - 55 2 ) 
- IsfSi (185 3 + 73S 2 3 ) - ^S!S 3 (1685 3 + 315 3 2 ) - ^S 2 (35 3 - 22S 2 3 ) + y 5 X 7 5 2 2 5 3 

i ^ C8^fiQ3 , q2\ 39 9 9 Ql0c 2 ^ Olio ^ cl2o i ^ cu\ 

+ — ^1 (,»^2 + ^3) ~ y^l^2^3 - J^l ^2 - g^l ^3 - J^l ^2 + — ^1 j ■ 

(25) 

More compactly, ordering it according to the power of S 3 , we have 

{•••} = ^{9S 2 (S 2 + 2S 2 ) 2 (S 2 - 2S 2 f - 245153 (S\ + 175 4 5 2 - 125 2 5 2 - 45 3 ) (S 2 - 2S 2 ) 

+ AS 2 (9Sf + 2825 4 5 2 - 205 2 5 2 - 85 3 ) (5 2 - 25 2 ) - 165! 5 3 (315? + 405 2 5 2 + 45 2 ) + 3205 2 5 4 

(26) 

Note although that the simplification in this example is not so significant as the previous 
one, which is because the power here is smaller. 



IV. J/ip HADRON PRODUCTION 

A. J/tp hadron production at leading order in v 2 

A well known result in heavy quarkonium physics is the J / ip hadron production in color- 



singlet model |19|-p2||, one out of six Feynman diagrams at leading order in a s is shown in 
Fig. [TJ Based on results in Sec. [TI], we will reproduce the behavior of the exact result by 
simple analysis. 

It is easy to find that denominator of the amplitude in Fig. [I] is proportional 
to (s — M 2 )(u — M 2 ), thus the summation of all six diagrams is proportional to 
[(s - M 2 )(t - M 2 )(u - M 2 )] In Feynman gauge, squaring the summed amplitude and 
summing/average over polarization, we get 



F( 3 S?) oc '"'"'^ . . (27) 

V 1 ; {{s - M 2 ){t - M 2 ){u - M 2 )] 2 



F 6 (M 2 ,s,t,u) 
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FIG. 1: One out of six Feynman diagrams for g + g — > cc[ 3 S[ ] + g, where all momentum are 
outgoing. The other five diagrams can be obtained by permutating the three gluons. 



Note that, an additional dimensionless factor (O j ^( 3 S^)} /M 3 and a flux factor l/2s are 
necessary to get the cross section of J/i/j production, but for our purpose it does not matter 
and we will neglect them. Before further discussion, we will prove that the behavior of 
this process to the differential cross section is oc M 4 jp\ in large limit, namely 
s, t, u ^> M 2 . Note that, to demonstrate it, we will choose gauge other than Feynman 
gauge. 

We begin with the amplitude in Fig. [I] (we label it with a "1" to denote that it is the 
first diagram. The other five diagrams will be labeled with 2, • • • , 6.) 



A 1 = Tr 



-- + mh a C- + m)(^gT a ^ 3 ] 



z(^3 + f + mj | 
(fc 3 + f ) 2 - m 2 



-igT^Y 1 ) 



+ m) 



:-h> 



y) 2 — m 2 



-igT a2 ^ 



-i2g 3 mTi [T a 'T a2 T a3 ] Tr 



7 "(l. + m ) 7 M3. 



f; 3 + f + m 



-r 



k,-P 



(28) 



where /i, (a^) is the spin index (color index) of the gluon with momentum hi, a is the spin 
index of the pair ccpSj 1 ] and m = M/2 is the mass of heavy quark. To get the above result, 
P 2 = 4m 2 and P a = 1 have been used. Notice that, if one sets m = 0, trace of the gamma 



chain in the last line in Eq. (J28Q vanishes because there are odd number of gamma matrixes 
then, which implies Ai oc M 2 . Therefore, squaring the summation of all six amplitudes and 
contracting with summations of polarizations, we will get oc M 4 /p|. if all denominators 
do not vanish when M — y 0. The only possible vanishing denominator comes from the 
summation of spin polarization of ccpS^ 1 ']: —g aa ' + P 1^T , where the second term violates 



This is possible because we will contract A\ with the summation of polariztion 



n aa' , rr 
9 + TF 
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the above argument. Thus, we should prove that contracting with M2 will at most give 
contributions at 0(M 4 ). 

Contracting the trace of the gamma chain in the last line in Eq. fl28|) with S , we get 



2 Tr 



+ Tr + m „ -^ 9 -^+m 



M3 ^3 + 



^2 



-7 



Ml 



hx-P 



-7 



M2 



(29) 



+ 0{M 2 ) . 



2k 3 ■ P 1 2k 2 ■ P 

It is convenient to choose a axial gauge so that P ■ e(fcj) = (z = 1, 2, 3), where e{ki) is the 
polarization vector of gluon with momentum hi. This choice of gauge does not change our 
argument in the last paragraph because it does not introduce small denominators. In this 
gauge, terms more than two f in Eq. (|29|) contribute to higher order in M 2 , thus 



At 



-2 Tr 



fJ-3 _ 



Tr 



Tr 



2k 3 ■ P 1 2k 2 ■ P 

$2 



+ 0(M 2 



2k 2 -P 



7 



t'2 



+ Tr 



7 



/i3 . 



t 



2k, P 



3 7 Miy*2 



+ 0(M 2 



(30) 



$2 



U 



+ Tr 



. $3 



S 



+ 0(M 2 ) . 



By circulating the three gluons, we get two other amplitudes which are proportional to 
Tr [T ai T a2 T a3 ], with 

r t i r u i 

(31a) 
(31b) 

Hence, using relation {/(fci), ^} = 0, we have Ai + A 2 + A 3 = 0(M 2 ). Considering that the 
summation of the other three amplitudes has the same behavior, we find contracting with 
will give contributions at 0(M 6 ). Finally, we complete the proof that J^- oc M 4 jp\. 



A 2 


= Tr 


s 


+ Tr 


7 Mi^i 7 M2 7 ^3 


+ 0(M 2 


A 3 


= Tr 




+ Tr 


u 


+ 0{M 2 



pa pa 

M 2 



As ^ oc M 4 /p| at large p T limit, F( 3 sj 1] ) and F 6 (M 2 ,s,t,u) in Eq. fl27|) must be 
proportional to M 4 . Using results in Sec. [TI], we get the general form 

F 6 (M 2 , s, t, u) = M 4 F 4 (M 2 , s, £, u) = S?(so5| + x^Ss + x 2 SlS 2 + x A S{) , (32) 

where Xj (i = 0, 1, 2, 4) are dimensionless numbers. The above simple analysis indeed repro- 
duced the behavior of exact result (e.g. one can find it in [17|), which expressed in our form 
is 



F( 3 S?) oc 



02/02 OO o2 1 1 o4A 

[(s - M 2 )(t - M 2 )(u - M 2 )] 2 



(33) 
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B. Relativistic corrections for J/ip hadron production 



A more interesting application of results in Sec. || is to understand the relativistic cor- 
rection behavior of J/ip hadron production in large p? limit. We will study in the limit 
s, t, m> M 2 , then result at leading order in v 2 as shown in Sec. |1V A| is 



C2 C2 

^i 1] )oc^?. (34) 

We will demonstrate that, in large pt limit, relativistic correction does not change this 
behavior, therefore, relativistic correction term is proportional to the leading term. The 
proof includes two steps. Firstly, we show that the denominator of summed amplitude 
is not changed by relativistic correction. Secondly, we show that the large behavior 
Jp- oc M A /p T is also not changed by relativistic correction. Combining these two points 

and using the general form in Sec. |T|, it is straightforward to find that relativistic correction 

s 2 s 2 

term is proportional to -^t-. 

Before two steps arguments, we briefly introduce the relativistic correction. In Sec. |IV A| , 
we use the color-singlet model to calculate the J / ip hadron production, where a charm (anti- 
) quark pair with definite quantum number (ccpsj 1 ']) is produced in hard collision. Further 
more, we did the non-relativistic approximation that relative momentum between cc pair is 
zero there. The relativistic corrections in this paper refer to corrections by expanding the 
relative momentum to higher power. For definiteness, we denote the momentum of c and c 
as 

' Pc = f + '' (35) 
Pc = f - q , 

where P is the total momentum of cc pair and q is half of the relative momentum. As 
the J/ip is a S'-wave quarkonium, only terms with even power of relative momentum in 
amplitude level contribute. In additional, projecting amplitude to S'-wave is equivalent to 
contract coefficients of relative momentum with terms like — g^ u + P ] J > 2 . 

Now, we are ready for the first step. The only possible source that may change the denom- 
inator of summed amplitude is the expansion of relative momentum for denominator with 
finite relative momentum. As an example, we study the expansion of upper denominator in 
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Fig.0 



1 1 _ 2k£ 



When contracting k% with —g^ u + pi ^K to get S-wave, p ^ gives leading contribution while 
— is suppressed as M — > 0. Notice that contracting fcg with ^ cancels the denominator 
(k 3 -P)~ 2 by one power exactly, therefore, we find that expanding this denominator to higher 
order in q does not change the denominator at large pt limit. This argument can be easily 
extended to any denominator that expanding to any power of q. Thus the first step is 
achieved. 

Let's then argue that the large pt behavior oc M 4 jp\ holds to all order of relativistic 
corrections, that is, terms of order M 2 jp\ vanish at the cross section level. For a finite 
relative momentum, the amplitude in Fig. |l| is 



Ax =Tr 

with 



{~\ + + ™)l a (^ + i + m)U(P, q, m) 



(37) 



U(P,q,m) = H,T 3 ^) {h+ P + q)2 _ m2 (-wT V 1 ) ^ _ P _ - g)2 _ m2 HgT V 2 ) • 

(38) 

As we study the large limit, both P and g have been boosted to a similar direction, say 
P, thus they have the following decomposition 

(39) 

q » = Lp» + + X 2q -^±n» , 

where £ = ^ = n is a light like momentum which satisfies (P ■ n) 2 /(n ) 2 ^> P 2 , and 
A is used to denote the power counting of corresponding term, that is, term proportional to 
A* behaviors as 0(M l ). Using the on-shell relations (±y + q) 2 = m 2 and the fact that 



(±- + + m)t(±j + ^ + m) = (±1 + C)^ + (±y + + m) , (40) 



we can rewrite Ax as 

1 (l - c 2 )P +2 



(— 2 + i + m )# 1 (-^ + J# + m h a (^ + </L + m)lft( — + j/ + m)il(P, g, m) 

(41) 
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where we also inserted two unit matrixes. Doing the Fierz transformation 



ha 



ij^-kl 



7 E T ii T ^ k i 



(42) 



with T A = 1,7 s , 7", 7V, ^7^2 and T x = T At , Ai becomes 



Ax 



4(1-C 2 )P 



- VTr 

■2\ p+2 



rA(-^- + ^ + m) 7 Q (- + ^ + m) 



x Tr 



- V Tr 



2(1 _ C , )P « a /_ 2 s 

x Tr 



p p 
f A (-| + ^ + m) 7 Q (| + ^ + m) 

P P 

+ + m)f A (^- + <| + m)n(P, g, m) 



where relations 



T A <g> ifiT x r/i = , , 2^ <g) ^ , 27 5 ^ ® 7^ , 27^ <g> 7^ 



(43) 



have been used, and f 1 



rfifi , f 3 = 7^. Note that we can think 



Tr 



f A (— y + ^ + m)7 a (-^ + ^ + m) to be O(M) effectively because squaring it will give 



a 0(M 2 ) result, thus we do leading power approximation for other terms and gives 



gp+2 



A=l,2,3 



T X {-- + i + m)l a (- + i + m) 



Tr 



^fYn(A^Ao) 



0(M 2 ) , 
(44) 



with 



i+C 



n(P,-A0) = ( —igT a ' A 7^ 3 ) — ^— ( —igT ai 7^ — — 



A: 2 -¥P) 2 ' 



-igT a2 ^ 2 ) . (45) 



Similar as Sec. IV A , we choose a axial gauge so that P ■ e(ki) = 0, then some terms in 
I1(P, |P, 0) do not contribute because there is a f on its either side. Thus, 



n(p, ^p,o)~(-z^T a v 3 ; 



*0 



i-C 



-igT ai ^ 



-h - ^py 



-igT a2 Y 2 ) 



(46) 



"1 7 M. 



(Jfe 3 + P) 2 (fc 2 + ^) 2 



2 
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and 



A 1 



2(1 



C 2 )P+ 2 

x Tr 



Tr [T ai T a2 T as ] ^ ^ Tr 

A=l,2,3 



- + 4 + m) 1 a (- + 4 + m) 



7 



-$ 2 



(A* + P) 2 (fca + P) 



:7 



A' 2 



(47) 



+ 0(M 2 



Observing that A = 3 does not contribute because there are odd number of Dirac matrixes 



in the last trace; A = 2 does not contribute because Tr 
e nPqa fay^ ^4 i s odd in in q. As a result, we get 

,3 



f 2 (-| + 4 + m) 7 a (f + d + m) 



oc 



-iff- 



2(1 -C 2 )P+ 2 

x Tr 



Tr p^y^y^J Tr 



^(-l7 + ^ + ^)7 Q (^ + ^ + ^) 



7 



MS . 



/'I . 



(h + py 7 (k 2 + p) 2 



7 



2(1-C 2 )P+ 



Tr p -1 " 1 T" 12 T a3 ] Tr 



__ + + m ) 7 < 



x <^ Tr 



n $2 



+ Tr 



s 



+ 0(M 2 ) 



0(M 2 



(48) 



Comparing it with Eq. (|3~0D and using similar argument, we find Ax + A% + A3 = 0(M 2 ) 
and that ^ oc M 4 /p| holds. 

Combining the above two points, we know that, in large px limit, the denominator is S^ 2 
and the numerator has a factor Sf. The rest numerator, as it is symmetric in s, t and u and 

has mass dimension of [M] 8 , must be proportional to 5*1 using the general form in Sec. ||. 

s 2 s 2 

Finally, we find that relativistic correction term is proportional to which is the same as 
term at LO in v 2 . This behavior for NLO relativistic correction has been found in Ref. [IT) , 
but the generalization to all order in this work is new. 

We conclude this section by explaining the logic to prove oc M A jp\. After the Fierz 
transformation, we in fact factorize the amplitude to soft parts (with scales of O(M)) and 
hard parts (with scales of 0(pt))- This factorized form is equivalent to double parton 
fragmentation formula in Ref. ||, which give contribution of 0(M 2 ) at cross section level. 
Factorized terms are then shown to vanish at large pt limit. Therefore, remained terms can 
only give contribution at 0(M 4 ). 
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V. SUMMARY AND OUTLOOK 



In heavy quarkonium physics, many processes involve three gluons, which result in s, t, 
u symmetric cross sections or decay widths. In this work we study general form of s, t, u 
symmetric polynomials, and find that they can be expressed as polynomials of Si, S 2 and 
S3 where the symmetry is manifest. For massless case the general form is summarized in 
Eq. (|T8|), and for massive case the general form is summarized in Eq. (|I§D . These general 
forms can be used to simplify expressions that are symmetric in s, t and u. Asymptotically, 
this method can reduce the length of original expression to one-sixth. Based on these 
general forms, one can also predict many interesting results by simple analysis. We give two 
examples regarding J / ip hadron production in this work. In the first example we work within 
the color-singlet model at LO in v 2 . By only arguing that the differential cross section has the 
behavior J^- oc jp\ in large px limit, namely s, t, u S> M 2 , we successfully reproduce 
the exact differential cross section up to four unknown constant numbers. In the second 
example, we consider relativistic corrections for color-singlet model in large pr limit. By 
showing that, in large pt limit, relativistic corrections do not change denominator and the 
behavior oc M jp%, we prove that differential cross section proportional to holds to 
all order in v 2 . This proof not only explains the proportion relation at NLO in v 2 found in 



Ref. [IT], but also generalizes it to all order. In view of that the s, t, u symmetry can give 
so many constraints, it is also interesting to study symmetry induced by four or even more 
gluons. 

Proportion relations at NLO in v 2 are also found for color-octet channel [IB]. However, 
since differential cross sections for color-octet channel have the behavior Jp- oc M 2 jp\ or 
Jp- cx 1/pt, proportion relations for them can not be constrained by only s, t, u symmetry. 
We will study this problem in a forthcoming work [[23]. 
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The Feynman diagrams were drawn using Jaxodraw 
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